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Abstract 

Singularities in the thermodynamics of Kerr-Newman black holes 
are commonly associated with phase transitions. However, such in- 
terpretations are complicated by a lack of stability and, more sig- 
nificantly, by a lack of conclusive insight from microscopic models. 
Here, I focus on the later problem. I use the thermodynamic Rie- 
mannian curvature scalar R as a try to get microscopic information 
from the known thermodynamics. The hope is that this could facil- 
itate matching black hole thermodynamics to known models of sta- 
tistical mechanics. For the Kerr-Newman black hole, the sign of R is 
mostly positive, in contrast to that for ordinary thermodynamic mod- 
els, where R is mostly negative. Cases with negative R include most of 
the simple critical point models. An exception is the Fermi gas, which 
has positive R. I demonstrate several exact correspondences between 
the two-dimensional Fermi gas and the extremal Kerr-Newman black 
hole. Away from the extremal case, R diverges to -|-oo along curves 
of diverging heat capacities Cj^# and Cn,Q, but not along the Davies 
curve of diverging Cj^q. Finding statistical mechanical models with 



* rupp einer @ ncf . edu 



like behavior might yield additional insight into the microscopic prop- 
erties of black holes. I also discuss a possible physical interpretation 
of \R\. 

Suggested PACS Numbers: 04.60.-m, 04.70.Dy, 05.40.-a 

1 INTRODUCTION 

A Kerr-Newman black hole is characterized solely by its mass M, angular 
momentum J, and charge Q [T]. Such simplicity allows a thermodynamic 
representation with laws analogous to the standard laws of thermodynamics 
[21 El m [5]. Previously, I discussed this structure in the context of ther- 
modynamic fluctuation theory [6]. This leads naturally to thermodynamic 
Riemannian geometry; see [7] for a review. 

The resulting thermodynamic Riemannian curvature scalar R has been 
explored by a number of authors for black hole thermodynamics |B1 El [IDl [HI 
[nillSIllIlSlIIHlinillHlITniEniiniEaE^ The main contribution 

of the present paper is an attempt at a physical interpretation of R, and its 
systematic evaluation for the Kerr-Newman black hole. The analogy with 
ordinary thermodynamics is emphasized, as is the significance of the sign of 
R. 

The thermodynamic fluctuation formalism requires stability, namely, fluc- 
tuations about a maximum in the total entropy. This issue poses difficulty 
for black holes. In ^ stability was obtained by restricting the number of 
independent fluctuating variables. In addition, an infinite extensive envi- 
ronment was employed to have the fiuctuations depend only on the known 
thermodynamics of the black hole. 

For an ordinary thermodynamic system, \R\ was interpreted p6] as pro- 
portional to the correlation volume where d is the system's spatial dimen- 
sionality and ^ is its correlation length. Direct calculations in a number of 
statistical mechanical models have verified this ; see [TT] for a more recent 
review. A thermodynamic quantity, R, then reveals information normally 
thought to reside in the microscopic regime, ^. Thus, R has been of interest 
in black hole physics, which has good thermodynamic structures, but little 
conclusive microscopic information. 

I interpret |-R| for black holes as the average number of correlated Planck 
areas on the event horizon. Although I give no direct microscopic model 
evidence, this interpretation would seem to be well motivated by the analogy 
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on the event horizon fluctuating independently of each other. Diverging 
\R\, which I take as signalling a phase transition, indicates highly correlated 
pixels. 

R diverges for extremal Kerr-Newman black holes, where the temperature 
T ^ 0. I demonstrate here, and previously [21], several exact limiting results 
matching extremal Kerr-Newman black hole thermodynamics to the two- 
dimensional (2D) Fermi gas {d = 2). Two dimensions are consistent with the 
membrane paradigm of black holes [30J. 




I also find instances of diverging R where the heat capacities Cj$ and 
Cq^q diverge and change sign, signalhng a change of stability. Although 
divergences in these heat capacities were identified by Tranah and Landsberg 
|31j in 1980, they have been little discussed in the literature. In contrast, I 
find no diverging R along the Davies curve [4j where the more familiar heat 
capacity Cj^q diverges. 

This paper is organized as follows. First, I review the thermodynamic 
fluctuation picture in [6]. Second, I discuss thermodynamic Riemannian 
geometry and curvature, including an attempt at a physical interpretation 
for R. Third, I calculate R for the Kerr-Newman black hole. Fourth, I 
compare the results with those in ordinary thermodynamics. 



2 THERMODYNAMIC FLUCTUATION THE- 



A major element in my approach is that the black hole resides in an infinite 
environment, characterized by mass Me, angular momentum Je, and charge 
Qe- The thermodynamics of the environment should be extensive; namely. 
Me, Je, and Qe should each scale up in proportion to the environment's 
volume. With this structure, thermodynamic fluctuations require only the 
known thermodynamics of the black hole. The environment's thermody- 
namic properties, which might be difficult to determine (dark matter, etc.), 
only sets the state about which fluctuations occur. 

^Such an interpretation is consistent with the assumption "that all the statistical de- 
grees of freedom of a black hole live on the black hole event horizon" [27] . 

^See [28] for a semipopular discussion picturing the quantization of area on the event 
horizon in terms of Planck areas. 
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This structure is thermodynamically unstable if we allow an exchange of 
all three variables (M, J, Q) [6l |3T]. Stability requires either a finite environ- 
ment or a restriction on the number of independent fluctuating variables. In 
[6] I took the later approach, and considered the stability of seven cases in 
an infinite environment: fluctuating (M, J, Q), (J, Q), (M, Q), (M, J), M, J, 
and Physically, we imagine that one (or two) of M, J, or Q is so slow 
fluctuating that we can consider it to be essentially fixed. 

I use geometrized units with M and Q in cm, and J and entropy S in 
cm^ PJ. Useful are the Planck length 



Lp = = 1.616 X lO^^^cm, (1) 



and the Planck mass 

/ lie 

Mp = y-= 2.177 xlO-^g, (2) 

with h, c, and G the usual physical constants. In geometrized units G = c = 1 
and Lp = Mp. 

The entropy of the Kerr- Newman black hole is ^ [33] 

1 



S{M, J, Q) = g [2M^ -Q^ + 2^M4 - J2 _ m^Q^J . (3) 
To convert S to real units, where it is Sbh, use 

with ks Boltzmann's constant [6]. 
The total entropy of the universe is 

Slot = Sbh + Se, (5) 

•^Kaburaki et al. '32] considered the stability of the Kerr black hole {Q = 0) under 
a variety of conditions, including ones similar to those here. These authors used the 
sophisticated Poincare turning point method which allows stability statements for finite 
environments with thermodynamics not necessarily known. With an infinite environment, 
however, the black hole entropy Hessian determinants are sufficient for considerations of 
stability [6]. Fluctuating conserved variables M, J, or Q are commonly referred to as 
"canonical," and fixed conserved variables as "microcanonical." 
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where is the entropy of the black hole's environment. The fluctuation 
probability is given by Einstein's formula [M] , 



Poce.p(|i). (6) 

Introduce the notation 

{X\X^X')^iM,J,Q) (7) 

and 

Fa = (8) 

with corresponding properties of the environment denoted by the subscript 
e. The intensive Fg^ values are independent of the size of the environment. 

Let us assume (incorrectly, as it turns out) that the black hole and the 
environment are fully in equilibrium, with a local maximum for Stot- Consider 
a small fluctuation AX" away from this equilibrium. Expanding each of the 
entropies in Eq. ([s]) to second order yields 



1 dF IdF 
AStot = F,AX^ + F,,AX: + -^^AX^AX^ + -^AX^AX:, 



(9) 



where the coefficients are evaluated at the equilibrium state, which is set by 
the environment. The conservation laws demand 

AX" = -AX:, (10) 

and a necessary condition for maximum entropy is 

Fa. = Fea. (ll) 

With a very large environment, the second quadratic term in Eq. (|9| is 
negligible compared with the first. To see this, fix the values of AX°, which 
equal —AX". As the environment is scaled up to infinite size at fixed F^a, 
X" scales up in proportion without limit, and dF^a/dX^ —>■ 0. The ability 
to drop this second quadratic term is a significant simplification offered by 
an infinite, extensive environment. 
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Eq. ([9]) now can be written as 
AStot 1 



--^?^,AX^AX^ (12) 
Kb ^ 

where the symmetric matrix^ 

^^^ = ~[l^J dX-dXP- ^^^^ 
The Gaussian approximation to the fluctuation probabihty is 

PdX^dX^dX"" = exp (^-^g^^AX^'AX"^ dX^dX^dX"", (14) 

where If?] is the determinant of Qa/s and n = 3 is the number of independent 
fluctuating variables. If we set one or two AX°'s to zero, reducing the value 



of n, Eqs. (12) and (14) are only trivially modified. Entropy maximum 
requires that the matrix of the remaining variable(s) be positive definite. 
Complete discussion of this is given in [B] . 



The first fluctuation moments are zero [31]: 

(AX°) = 0. (15) 
The second fluctuation moments are 

(AX°AX^) = c/"^, (16) 



with (yf"^ the components of the inverse ga/s matrix. 

Further notation is given in the Appendix, where I define the basic ther- 
modynamic variables T, Q, and $, the simplifying variables a, (3, K, and L, 
the entropy Hessian determinants p2, P2, and pg, with numerators A, B, and 
C, and the heat capacities Cj,q, Cj,$, and Cn,Q. Diverging heat capacities 
are important below, and Figure 1 shows curves of infinities as well as the 
extremal limiting curve where the temperature T — > 0. 

"^In the symbol ja(3 was used for this quantity, ga/3 denoted this quantity without 
the unit conversion factor for S. 
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3 THERMODYNAMIC RIEMANNIAN GE- 
OMETRY 



In this section, I summarize the thermodynamic Riemannian geometry. 
3.1 Thermodynamic metric 

The metric follows naturally from the observation that the quadratic form in 



Eq. (12) transforms as a scalar under any coordinate change because AStot 



depends only on the initial and final thermodynamic states. Hence, 
9A 

(A/)2 = = g^.AX'^AX'' (17) 

Kb 

is a Riemannian line element. It is unitless and positive definite assuming 



stability. Its physical interpretation is clear from Eq. (14): the less probable 



a fluctuation between two states, the further apart they are. 



In the definition of ga/s in Eq. (13), S was converted to Shh/ks in real 
units, essential in Eq. ([6]). Such a unit conversion is unnecessary if R is not 
needed beyond a proportionality constant. However, a quantitative interpre- 
tation of R in analogy with ordinary thermodynamics requires a unitless line 



element of the form in the exponential of Eq. (14). 



To get the metric elements in Eq. (13), I used the special properties of the 



conserved variables (M, J, Q). Once we have Eq. ([13]), gap transforms as a 
second rank tensor under a change of coordinates j7]. Generally, under such 



a transformation the Hessian form in Eq. (13) will not persist. However, 
since we know the function S = S{M, J, Q), there is no need in this paper to 
change coordinates. 

3.2 Thermodynamic curvature 

Calculate R as follows [35] : the Christoffel symbols are 

r^T- = ^g"" {g^Pa + 9^.^,p - gp^,f.) , (is) 

where the comma notation indicates partial differentiation. The Riemannian 
curvature tensor is 

pa -pa -pa , -pfJ. -pa -pa fTO\ 
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and the Riemannian curvature scalar is 



R = g^'Rl^,. (20) 

R is independent of the choice of coordinate system, suggesting it is a fun- 
damental measure of thermodynamic properties. Since the line element is 
unitless, R will be unitless. 

For two-dimensional thermodynamic geometries (n = 2), all components 
of the Riemannian curvature tensor may be expressed in terms of the cur- 
vature scalar R [35]. Not so in higher dimensions. However, it was argued 
[36] that R is the essential quantity in thermodynamic geometry regardless 
of the number of independent thermodynamic variables. 

For an ordinary pure fluid, a common picture [26j is that of an open 
subsystem with fixed volume V surrounded by an infinite environment of the 
same fluid. The entropy fluctuation is 

kB ~ 2^ ksdx^^dx-^ ' ^ ' 

where s is the entropy per volume (in units of ks per volume), and x^ and 
are the energy and particle number per volume, respectively. The pure 
fluid line element was written Eni with V omitted: 



(A/)2 = -LJ^Ax'^Ax", (22) 
^ ^ ksOx'^dx" ^ ' 

and has units of inverse volume. The corresponding R has units of volume. 

Logically, however, the pure fluid line element could have been written in 
the unitless form 

with neither the subsystem entropy S nor the conserved energy and particle 
number {X^, X^} divided by the constant V . The form of this line element 



matches that of the black hole line element Eq. (17). It has unitless R. 



For the pure fluid, it was noted [26] that R calculated with the line ele- 



ment Eq. (22) is zero for the pure ideal gas, suggesting that i? is a measure 



of intermolecular interaction strength. Indeed, calculations showed \R\ to be 
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proportional to the correlation volume ^'^ for a number of statistical mechan- 
ical models [71 ITT] . 

Such calculations dovetailed nicely with R having units of volume. How- 
ever, the units of R are not naturally determined]^ With the equally valid 



line element in Eq. (23), the fixed V now appears in the denominator of R, 
and R is unitless. If we imagine the fluid broken up into three-dimensional 
{d = 3) pieces each of volume V, \R\ is the average number of correlated "pix- 
els." The physical interpretation of \R\ is then essentially the same regardless 
of whether or not we pull V out of the line element. 

This leads to a possibly useful way to look at black holes. Although there 
is no flxed subsystem volume to set a scale, the Planck length Lp suggests a 
physical constant for this rolej^ Black hole thermodynamics takes place on 
the 2D event horizon. It is natural to break it up into square pixels each of 
area [28]. By analogy with the pure fluid, I interpret \R\ as the average 



p 

number of correlated pixels. Figure 2 illustrates this physical interpretation. 

I cannot presently support this idea with microscopic calculations of a 
type which were so valuable in ordinary thermodynamics. However, the 
correspondence in the extremal limit with the 2D Fermi model [29], discussed 
below, indicates at least consistency with the black hole membrane paradigm 
[30] which puts all black hole properties on the 2D event horizon. 

Note, this interpretation is only of |-R|. Janyszek and Mrugala [3H] argued 
that the sign of R is also important. I amplify on this in Sec. 6. 

Finally, in a coordinate system with metric elements of Hessian form, R 
simplifles. The arguments in [36j allow one to show that with metric elements 



in Eq. (13) 



R = \g^''g^°g''' i9,u,^9on,p - g,o,pg.^,.) ■ (24) 

The second derivatives of the metric elements cancel in the calculation. 



^My previous arguments about the significance of volume units for R (see, e.g., Sec. 
VLB of |7]) may have been overstated. Model calculations and the path integral approach 
to thermodynamic fluctuation theory [S^ are the best way to establish R oc S,"^. However, 



the pulling out of V in Eq. (22 1, and the resulting units of volume for R, is certainly 
natural and leads to correct results. 

^Note, in ordinary thermodynamics there is no physical constant with units of length. 
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3.3 Background on black hole thermodynamic curva- 
ture from the entropy metric 

Aman et al. [16] presented a recent review of thermodynamic curvature in 
the context of black holes, so my remarks in this section will be brief. Ferrara 
et al. [8] were the first to apply thermodynamic curvature to black holes, to 
calculate critical behavior in moduli spaces. Cai and Cho [9] connected phase 
transitions in Bafiados-Teitelboim-Zanelli (BTZ) black holes to diverging R. 
They also identified a correspondence with R for the Takahashi gas, suggest- 
ing that an appropriate black hole statistical model might be a system of 
hard rods. 

Aman, Bengtsson, and Pidokrajt [TO] were the first to evaluate R for vari- 
ous instances of the Kerr- Newman black hole, especially the two-dimensional 
{n = 2) Kerr and Reissner-Nordstrom cases. They also considered a nonzero 
cosmological constant. Arcioni and Lozano-Tellechea [I2] worked out five- 
dimensional black holes and black rings, including an extensive review. They 
connected phase transitions to both diverging R and diverging second fluc- 
tuation moments. 

Aman et al. [T3| examined R in the context of homogeneous functions, 
emphasizing in particular cases with R = 0. Aman and Pidokrajt [I5] in- 
vestigated Kerr and Reissner-Nordstrom black holes in spacetime dimensions 
higher than four. They found that patterns in four dimensions continue to 
higher dimensions. Sarkar et al. evaluated R for a general class of BTZ 
black holes, including quantum corrections to the entropy. Mirza and Za- 
maninasab [18j worked out the curvature of the full 3D geometry for the 
Kerr-Newman black hole. They found that R diverges at the extremal limit, 
but not along the Davies curve. Aman et al. [20j reported results on dilaton 
black holes. 

3.4 Curvature from other than entropy metrics 

One is certainly not constrained to do thermodynamic Riemannian geometry 
with fluctuations and its entropy metric. Another possibility is to express 
the internal energy in terms of its natural variables, M = M{S, J,Q), and 
construct an energy metric from its Hessian. This was done originally in or- 
dinary thermodynamics by Weinhold [39]. The energy metric is conformally 
equivalent to the entropy metric [IQ] , with the same angles between vectors, 
but different R. 
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A motivation for exploring other metrics (including the energy metric) is 
a concern by some authors about the physical validity of cases with R = 
from the entropy metric. If R is interpreted as a measure of interactions 
among gravitating particles, one might logically expect \R\ to be uniformly 
large for black holes, where gravitational forces are very big. 

However, the interpretation of R in this paper takes a different approach. 
The gravitating particles have presumably collapsed to the central singular- 
ity, shrinking the interactions between them to zero volume. The statistics 
underlying the thermodynamics are envisioned to be on the event horizon. A 
result R = seems now not so unreasonable. Yet, I present little in the way 
of microscopic evidence for this point of view, so concerns about the physical 
validity of = certainly cannot be dismissed. 

For the Reissner-Nordstrom black hole Aman et al. [10] found i? = with 
the entropy metric. To avoid this zero curvature, Shen et al. [13] constructed 
a new entropy metric, replacing Q with $ and M with M — <^Q. These 
authors also argued that R should signal (by diverging) a phase transition 
at the Davies curve. Their modified R shows such a divergence. A detailed 
analogy with the van der Waals phase transition was worked out with their 
modified metric. The authors also connected to modern themes in particle 
theory, such as holography and the AdS/CFT correspondence. Mirza and 
Zamaninasab [IS] evaded zero curvature by evaluating R for the full 3D 
Riemannian geometry. Here, R is always positive, as will be discussed in 
Sec. 5.1. Quevedo and collaborators [19], EH ES] suggested that this issue 
requires Legendre invariant metrics to deal with properly. They constructed 
a detailed framework based on this idea. Medved [21] also gave a recent 
discussion of these issues. 

4 Background on phase transitions 

Exactly what constitutes a black hole phase transition is somewhat unsettled 
in the literature. In ordinary thermodynamics the modern belief is "a phase 
transition occurs when there is a singularity in the free energy or one of its 
derivatives" jH]. Phase transitions can bring about dramatic contrasts, like 
between a solid and a gas. Or changes can be more subtle, hke the onset of 
a gradual deformation in crystal structure. Conjectured phase transitions in 
Kerr-Newman black holes are typically of the more subtle variety, second- 
order phase transitions associated perhaps with a diverging heat capacity. 
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Phase transition theory in ordinary thermodynamics typically includes 
equilibrium between system and environment. Achieving this with black 
holes can be difficult. A more serious problem is the lack of conclusive mi- 
croscopic models for black holes. This makes it hard to identify objects as 
fundamental to phase transition theory as order parameters and correlation 
functions. 

There are then a number of viewpoints of what might be involved in a 
black hole phase transition: 1) a change in topology, 2) a divergence of a 
second fluctuation moment, 3) a divergence of a heat capacity, 4) an on- 
set of instability like that in an axisymmetric rotating self-gravitating fluid, 
5) a divergence of R, and 6) consistency with the scaling laws of critical 
phenomenal^ 

For the Kerr-Newman black hole, there is no topology change except 
possibly at the extremal limit. Otherwise, the topology is that of the sphere 

m- 

Second fluctuation moments are connected to quantities such as heat 
capacities through thermodynamic fluctuation theory [M], so viewpoints 2) 
and 3) are related, a point not always clear in the black hole literature. 

Davies [4j argued that the curve of diverging Cj^q constitutes a second- 
order phase transition. However, this has been disputed by a number of 
authors. One issue is whether or not the Davies curve marks an actual change 
of stability. I found that it does so only for M fluctuations [B]. Davies |1] 
also brought up the analogy with the change of symmetry of an axisymmetric 
rotating self-gravitating fluid [43j. However, this was questioned ^44j since the 
nonrotating black hole also crosses the Davies curve as charge is increased. 

Arguments based on scaling theory usually involve attempts to introduce 
an order parameter. Cai et al. [45j and Kaburaki [IH] argued that the 
extremal limit constitutes a second-order phase transition and proposed the 
difference between the inner and outer black hole radii as the order parameter. 
Lousto IITJ HH] emphasized instead a phase transition along the Davies curve. 
He used f2 — f2c as the order parameter, where VL^ is the angular velocity along 
the Davies curve. He also worked out critical exponents and discussed them 
in the context of scaling theory. Lau |49j also argued that the Davies curve 
corresponds to a second-order phase transition. 

^Another possibility is the first-order phase transition consisting of a black hole con- 
densing out of the background space However, this differs from the second-order phase 
transitions at issue here. 



12 



5 KERR-NEWMAN THERMODYNAMIC CUR- 



VATURE 

In this section, I work out R for (M, J, Q), (J, Q), (M, Q), and (M, J) fluc- 
tuations. M, J, and Q fluctuations, with n = 1, have i? trivially zero, and 
require no special consideration]^ 

I go beyond Aman et al. |10], and report all situations. In the Kerr- 
Newman family, these authors focused primarily on Reissner-Nordstrom black 
holes, represented by the geometry of (M, Q) fluctuations with J = 0, and 
Kerr black holes, represented by the geometry of (M, J) fluctuations with 

g = o. 



5.1 (M, J, Q) fluctuating 



Here, all {M,J,Q) fluctuate. By Eq. (17), 



2_ (7ii(AM)2 + 2(7i2AMAJ + 2(7i3AMAQ+ 
^"^'^ - g22iAjy + 2g,,AJAQ + g,,{AQf, ^^^^ 

corresponding to a 3D Riemannian geometry (n = 3). This case falls entirely 
outside the domain of stable fluctuations [6], |3T] , and so I give it only a little 
attention. 



Evaluation with Eq. (20) shows R to be always real and positive, with 
a minimum of {Mp/2^/^TMy at the origin J = Q = 0. R is shown in 
Figure 3|^ It has no anomalies except at the extremal limit, where it diverges 
proportional to T^^. 

Mirza and Zamaninasab [IH] also worked out this case. With zero cos- 
mological constant, they found that R diverges at the extremal limit, but 
nowhere else. In particular, they found no divergence along the Davies curve. 
They also found nonzero R for the Reissner-Nordstrom case, J = 0. Fig. 3 
corroborates these findings. 



^Trivial geometries {n ~ 1) reflect noninteracting systems. For example, open fluid 
systems characterized by one fluctuating parameter, usually the internal energy or the 
temperature, generally do not have interactions, e.g., a gas of photons. With interactions, 
an additional parameter, like the density, is generally required. 

^ As was pointed out [lOl [18] , the full algebraic expression for R is not particularly 
revealing. I will adhere to tradition and not show it here. 
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5.2 (J,Q) fluctuating, M fixed 



Here, {J,Q) fluctuate at fixed M. By Eq. (17), 

{Aiy = g22{AJf + 2(723AJAg + to(AQ)2, (26) 

corresponding to a 2D Riemannian geometry (n = 2). Fluctuations in this 
case are stable for all states in the physical regime [6]. 
By Eq. (20), and the definitions in the Appendix, 

{K^ + L'^K^ - 2K^ - 2K^ + 3L^K - 3K + 2) f Mp 



^ = ' [JT) ■ 

R is shown in Figure 4. As is argued in the Appendix, B is never zero in the 
physical regime. Hence, R only diverges at the extremal limit where K ^ 0. 

Let us examine further the extremal limit. Equations (27), (55), and (64) 
yield the extremal limiting expressions 

Cj,Q = ^M'L^T, (28) 

and 

2M2 

The limiting product of curvature and heat capacity, 

^ (^) (^) 

is a unitless, scale free constant independent of where we are on the extremal 
limiting curve. In Section 6.2, I evaluate the statistical mechanics of the 2D 
Fermi gas and demonstrate several exact correspondences with these results 
at low temperature. 

We have from Eq. (16), and the definitions in the Appendix, the dimen- 
sionless second fluctuation moments 

J- 9 ?r / u\ 

(31) 



((AJ)2) 1 Ik^ - L^K + 2K ( M 



h V27T^ B \Mp^ 
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and 



X(Ag)2) _ 1 /137.04 biK^ + L^K - K) 
e ~ 2V^^V B ' 



(32) 



where e is the electron charge and 137.04 = h/e^ is the fine structure con- 
stant. Both fiuctuation moments are real and nondiverging in the entire 
physical regime. They have maxima of l/v^27r and 3.302, respectively, at the 
origin J = Q = 0, and decrease to zero as -v/T at the extremal limiting curve. 



5.3 {M,Q) fluctuating, J flxed 



Here, (M, Q) fiuctuate at fixed J. By Eq. (17), 



(33) 



corresponding to a 2D Riemannian geometry {n = 2). There is a slice of 
stability [6] in the {y/a, y/]3) plane bounded by the extremal limiting curve 
and the curve C = 0, as shown in Figure 5. 



By Eq. (20), and the definitions in the Appendix, 



R 



(L-1)(L + 1) 



( 3KL^ - AL^ + AK^L^- 
8K^L^ - QKL^ + 36L^+ 
K^L^ - AK^L"^ + IAK^L'^+ 
AOK^L^ - 36KL'^ - 9QL^+ 
8K^ + AK^ - 3QK^- 
V 32K^ + ASK + 64 



\M ) 



(34) 



R is shown in Figure 6. Despite only a limited slice of stability, located in the 
"saddlebags" near J/M^ = ±1, i? is real everywhere in the physical regime. 

Along the line J = 0, we clearly have i? = 0, since L = 1. This was 
demonstrated by Aman et al. [TU] who also pointed out that there is no 
curvature anomaly at the Davies point {J/M'^^Q/M) = (0,0.8660). I add 
that no point along the line J = lies in the stable slice, as is clear in Fig. 
5. 

Figure 6 shows a steep drop to negative R near Q/M = ±1. The cut 
away view shows this as a waterfall shape. Such abrupt behavior, the only 
case of negative curvature for the Kerr- Newman black hole, reminds one of 
the abrupt change in sign for the Takahashi gas [50] and for the finite ID 
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Ising model [51], which will be discussed further in Section 6.1. None of these 
negative values fall into the stable slice, however. 

Closer to the extremal limit, R comes up again, diverging to +00 at all 



points on the extremal limiting curve K = except J = 0. Equations (34) 



and (55) yield the extremal limiting expression 



Despite the difference between (M, Q) and (J, Q) fluctuations, this limiting 



expression is the same as Eq. (29) for {J,Q) fluctuations, and the match to 
the 2D Fermi gas applies equally well here. 

R has an additional divergence, to +00, at the other boundary of stability, 



C = 0. R diverges as C ^. Cj^$ diverges as C ^, by Eq. (66) 



We have from Eq. (16), and the definitions in the Appendix, the dimen- 



sionless second fluctuation moments 



((AM)2) 1 ^M,\ lK^-Lm + 2K 

, (36) 



with rrif. the electron mass, and 



((AQ)2) 1 /137.04 -AK^ + 2L^K^ -8K^ + 2L^K , , 

. (37) 



e 2V vr V C 

Both these quantities are real in the slice of stability. They go to zero as ^JT 
at the extremal limit -ft^ = 0, and diverge along the curve C = 0. Hence, 
changing stability is marked by both diverging R and diverging fluctuations. 

Note, fluctuations in M expressed in units of the electron mass are huge. 
However, in units of the Planck mass they would be much smaller, on the 
order of the fluctuations in J and Q. 

5.4 (M, J) fluctuating, Q flxed 



Here, (M, J) fluctuate at fixed Q. By Eq. (17) 



(A/)^ = (7n(AM)^ + 2guAMAJ + (?22(AJ)^ (38) 
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corresponding to a 2D Riemannian geometry (n = 2). Stability is confined 
to a slice bounded by the extremal limiting curve and the A = curve, as 
shown in Figure 7. 



By Eq. (20), and the definitions in the Appendix, 



R 



1 



27rKA^ 



( + 3K^ + 2L'^K^ + QL'^K^- 

SL'^K^ + AL^K^ - 8K^ + QL'^K- 
V 21L'^K + 12K + 9L^ - 24^^ + 16 / 



M ) 



(39) 



It is shown in Figure 8. Despite only a limited slice of stability, R is 
real and positive everywhere in the physical regime. Its minimum value 
is {Mp/2^/7lMf at the origin. 

Aman et al. [10] computed i? for Q = 0, and found it to diverge at the 
extremal limit. They pointed out that there is no curvature anomaly at the 
Davies point (J/M^,Q/M) = (0.6813,0). This is confirmed by the findings 
here. I add that no point with Q = lies in the stable regime, as is clear in 
Fig. 7. 

R diverges at the extremal limit K = Equations (39) and (55) yield 
the extremal limiting expression 



R 



2Ml 



TiM^L'^T 



(40) 



Remarkably, this is the same as Eqs. (29) and (35) found previously. 



R has an additional divergence, to +oo, at the other boundary of stability, 
A = 0. R diverges as A~'^. Cq^q diverges as A~^, by Eq. (67). 



We have from Eq. (16), and the definitions in the Appendix, the dimen- 



sionless second fiuctuation moments 



and 



((AM)2) 



((AJ)^) 

h 




+ L'^K - K 



A 




-2K^ + L'^K^ - AK^ + L^K f M\ 
A [Wj 



(41) 



(42) 



Both these quantities are real in the slice of stability. They go to zero as 
-\/T at the extremal limit K = 0, and diverge along the curve A = 0. Hence, 
changing stability is marked by both diverging R and diverging fiuctuations. 
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6 DISCUSSION 

In this section, I review evaluations of R in ordinary thermodynamics, and 
compare with the Kerr-Newman black hole. 

6.1 Curvature in ordinary thermodynamic models 

Table I reviews signs and divergences of thermodynamic curvature in several 
ordinary thermodynamic models. Most table entries are simple models where 
R may be worked out in closed form. The tendency is negative R where 
attractive interactions dominate, and positive R where repulsive interactions 



dominate, Janyszek and Mrugala p8] emphasized the importance of the 
sign of R, and identified the quantum gasses, 3D Bose and 3D Fermi, as 
essential examples with opposite signs. 

The signs of R for the standard critical point models in Table I are all neg- 
ative, and have critical point divergences R — > — oo. This is quite unlike the 
Kerr-Newman black hole, with its predominantly positive R and divergences 
R +00. 

Table I shows a group of weakly interacting systems with "small" \R\. 
Small means on the order of the volume of an intermolecular spacing or less. 
I view such values of R as physically equivalent to zero, since the meaning 
of correlation volumes of this size is lost in the "noise" of thermodynam- 
ics breaking down as individual atoms and spins become visible. The ID 
antiferromagnetic Ising model [52| 153] is perhaps misplaced here, since its 
disaligning interactions might propagate a long way. However, for antifer- 
romagnets, the true ordering field is a staggered field, and not the constant 
field used for the calculations in Table I. A reassessment of this model in 
these terms might be called for. 

There are three cases in Table I having both positive and negative curva- 
tures. The ID g-state Potts model [HI has sign related to q. For q > 2, 
and nonzero field, there are significant regions of negative R at low tempera- 



^"^One must guard against the impression that there is a connection between thermo- 
dynamic stability in the sense here and the sign of R. It is tempting, for example, to 
envision the n = 2 thermodynamic Riemannian geometry of the type here as a 2D surface 
embedded in a 3D flat Euclidean space from which it inherits its metric. In such a con- 
struction, thermodynamic stability requires R for the 2D embedded surface to be negative. 
However, this picture is incorrect, as has been discussed in Sec. IV. G of Ref. [7 . There is 
no connection between thermodynamic stability and the sign of R. 
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ture. The 2D Potts model has the dimensionahty of the Kerr-Newman event 
horizon. Its R has not yet been evaluated, but perhaps its study could yield 
an appropriate critical line with positive R. 

The Takahashi gas [50] has the typical negative R in the gaslike phase 
where attractive interactions dominate, and small in the liquidlike phase 
where interactions are short-range. However, going from one phase to the 
other by changing the density at constant temperature, there is a pseu- 
dophase transition accompanied by a sharp positive curvature spike. Cai 
and Cho [9j connected this spike to a phase transition in the BTZ black hole. 

An abrupt change in sign of R is also present in the finite ID Ising ferro- 
magnet of spins [M]. This model has the typical negative R for large N, 
but a sharp rise to large positive values as N is decreased. Whether or not 
this result has relevance here is unclear. 

At the bottom of Table I there are the 3D Fermi gas |38] and the 3D 
Fermi paramagnet [61]. Both models have positive i?, diverging as T — 0. 
These results lead me now to take a closer look at Fermi gasses, particularly 
the 2D Fermi gas. 



6.2 Curvature for the 2D Fermi gas 

For the 3D Fermi gas at low T, R seems to diverge 
as T~ 



^ i] as T~^/^, and not 

^ in Eq. (29) for the Kerr-Newman black hole. This motivates me to 

By the reasoning leading to Eq. (8.1.3) of [62] . 



work out the 2D Fermi gas 
the 2D Fermi gas has thermodynamic potential 



'T' 



T 



kB9>^ f2{v) 



(43) 



with pressure p, r] = exp{fi/kBT), chemical potential /i, thermal wavelength 
A = h/ y/2TTmkBT, particle mass m, weight factor g = {2s + 1), particle spin 
s, and 



1 



X dx 



(44) 



I use obvious fluid units for all quantities, including S and T. The integral in 
Eq. (44) converges for /2(?7), and yields /i(?7) = ln(l + rj). fo{v) f-i{v) 
follow from fi{ri) using the recurrence relation fi-i{i]) = TjfKrj) [62] . 
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Define the heat capacity at constant particle number and constant area 
A by 



Cn,a = T 




Nk 



B 



'fM foiv) 



(45) 



The second equality is by Problem S.lO.ii of [62j. The methods of [62] now 
yield the limiting low T expression 



C 



AkB 3/^2 
Evaluating R with Eq. (6.31) of [7] yields 



(46) 



R 



(47) 



Numerical evaluation over the physical range —oo < fi < +oo and < T < 
oo indicates R is always positive. The methods of [62j yield the limiting low 
T expression: 



R 



3/^2 



(48) 



2TT^gmkBT 

The limiting T dependences of Cn,a and R match the corresponding 



Kerr-Newman black hole quantities Eqs. (28) and (29). This connection to 



a 2D model is consistent with the membrane paradigm of black holes 
Furthermore, the limiting product of curvature and heat capacity, 



R\ f Cn^a 
a) [~k^ 



hi 



3h^ 



\2t:^ gmk bT A ^ 



'2TT^gmkBTA^ 
3/^2 



1, 



(49) 



is a unitless, scale free constant independent of density. The factor A below 
R undoes the traditional pulling out of A in the ordinary thermodynamic 
line element. R/A here is analogous to R for the Kerr-Newman black hole. 



The constant products Eqs. (30) and (49) are equal, remarkable for systems 
apparently so different. 

However, note a key difference. The Kerr-Newman black hole entropy 
Eq. ([3]) does not go to zero in the extremal limit, as it does for the 2D 
Fermi gas with its unique ground state. Resolution probably requires a more 
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sophisticated Fermi gas model. Note as well that I have presented no detailed 
correspondence between the Kerr-Newman black hole thermodynamics and 
a specific microscopic Fermi model. Such a connection is necessary to make 
the results given here something more than a possibly useful direction to 
explore. 

6.3 Curvature for black hole critical points at nonzero 

T 

For the phase transitions found at the non-extremal boundaries no appro- 
priate models with evaluated i?'s present themselves. The signs of R of the 
simple critical point models in Table I are all negative, in contrast to the 
Kcrr-Ncwman black hole results. Hence, I make no attempt here to suggest 
an order parameter or to evaluate and interpret possible critical exponents 
and scaling relations between them. 

7 CONCLUSIONS 

In conclusion, the following were done in this paper for thermodynamic Rie- 
mannian geometry based on the entropy metric. 

First, 1 attempted a physical interpretation of R for black holes. It was 
based on analogy with the interpretation in ordinary thermodynamics. Per- 
haps, this interpretation lessens concern over the physical plausibility of the 
occasional result R — 0. 

Second, I reviewed previous evaluations of R and phase transitions in 
Kerr-Newman black holes. 

Third, I gave a complete evaluation of R for the Kerr-Newman black 
hole. In all cases, R was found to be positive in stable fluctuation regimes 
and to diverge to -|-oo at the extremal limit. I also found it! to diverge to 
-|-oo at nonzero temperatures along curves of changing stability, where the 
heat capacities Cj,$ and Cq^q diverge. 

Fourth, I argued that the sign of R is important, and tabulated signs in a 
number of ordinary thermodynamic models. I found that most of the simple 
critical point models have negative R. This might make them problematic for 
understanding Kerr-Newman black hole phase transitions. Different models 
might be required. 



21 



Fifth, I noted that the Fermi gas is one of the few known cases in or- 
dinary thermodynamics with large positive R. I estabhshed several exact 
correspondences between the 2D Fermi gas and the extremal Kerr-Newman 
black hole. This suggests that microscopic models with fermions might be 
useful as a framework for formulating a microscopic description of black holes. 



Notation was defined in [6j, and is summarized here. I differ only with the 
metric elements Qa/s, including here the unit conversion factor for S in Eq. 



Define the temperature T, the angular velocity f2, and the electric poten- 
tial $ by H ED 



8 APPENDIX: NOTATION 




(50) 




(51) 



and 




(52) 



Two standard unitless variables are |1] 



{a,/5}^{JVM^gVM2}. 



(53) 



Simplifying the notation are [21] 




(54) 



We may show that 



1 _{K^ + 2K + L2) M 
f ~ AK 



(55) 
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To be in the physical regime of real S and T requires 



a + /5 < 1. 



(56) 



The curve of equahty, a + [3 = 1, has K = T = and constitutes the 
extremal limit, thought to be unattainable by the third law of black hole 
thermodynamics [63] . 

Major components in the discussion of stability are the entropy Hessian 
determinants^ 



P2 



P'2 



ZJP. I 



L 



2\ 2 



and 



L 



2\ 2 



9n 


912 


921 


922 


922 


923 


932 


933 


9u 


913 


931 


933 



-2K^ - 3K^ - 2L^K + 2K - + A 
+ L^K -K + 1 



2L'^K -AK - 2L^+ 
V 10L2 - 8 



,(57) 
(58) 



(59) 



The numerators of p2, P2, and P2 are, respectively, 
A = -2K'^ - 37^2 _ 2L'^K + 2K - 31"^ + 4, 
B = + L'^K - K + 1, 

and 



C 



L^K + 21"^ K -AK ~ 2L^ + lOL^ - 8 



(60) 
(61) 

(62) 



Curves along which these numerators go to zero identify changes of sta- 
bility accompanied by divergences of heat capacities. A = in the physical 
regime if and only if 



a 



(3-4/3)/j^ 
4(/3-l)2 



(63) 



^^In 6 the metric elements gai3 did not include the conversion factor for S, we must 
undo these to make the entropy Hessian determinants the same as in [6]. 
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This curve is shown in Fig. 7. B is never zero in the physical regime, since 
i^' > and L > 1. C = along a single curve in the physical regime, shown 
in Fig. 5, with its algebraic expression too complicated to show here. 
Finally, the heat capacities j31j 



_^(dS\ _ M^KjK^ + L^ + 2K) 
Cj,Q = ^\dT)j^- 4(L2 - 2K) ' ^^^^ 

(05) 

which evaluates to 

_ M^K [K^ + 2K + L^) ( -L' - K^L^ + KL^ + AL^+ \ 

AC \K^ + AK^ + 2K-2 J' ^^^^ 

ancS 

Cj^Q diverges if = 2K. This may be written 

+ 6a + 4/3 = 3, (68) 
which gives the Davies curve, shown in Fig. 1. 
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3D Fermi gas ^38j 
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Table I. Thermodynamic curvature for ordinary thermodynamic systems. 
I give the number of independent thermodynamic parameters n, spatial di- 
mension (i, sign of i?, and comment on possible divergences. In some systems 
d is not set, and I denote this with All signs of R have here been put 

into the sign convention of Weinberg [35]. An indication small" means 
\R\ has a value on the order of the volume of an intermolecular spacing or 
less. 
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Figure Captions 

FIG. 1. Some characteristic curves for the Kerr-Newman black hole; see 
the Appendix. The curve along which Cj^q diverges is the Davies curve. R 
diverges at the extremal limit and along curves corresponding to a change of 
stability, which have diverging Cj_$ and Cq^q. 

FIG. 2. The event horizon broken up into Planck area pixels. The dark 
pixels are portrayed as somehow correlated. I propose that \R\ measures the 
average number of correlated pixels. 

FIG. 3. R{M/Mpf as a function of J/M^ and Q/M for (M, J, Q) fluctua- 
tions. R is real, positive, and regular in the physical regime, and diverges as 

at the extremal limit. 
FIG. 4. R{M/Mpf as a function of J/M^ and Q/M for (J,Q) fluctuations. 
R is real, positive, and regular in the physical regime, and diverges as 
at the extremal limit. 

FIG. 5. Stable fluctuation regime for (M, Q) fluctuating at fixed J is in- 
dicated by -I- signs. The case with J — Q corresponds to the Reissner- 
Nordstrom black hole, which lies entirely out of the stable regime. 
FIG. 6. R{M/Mpf as a function of J/M^ and Q/M for (M, Q) fluctuations. 
R is real everywhere in the physical regime. It is mostly positive, but there 
are two regimes of negative values near Q/M — ±1. R diverges at both 
limits of stability. 

FIG. 7. Stable fluctuation regime for (M, J) fluctuating at fixed Q is indi- 
cated by + signs. The case with Q = corresponds to the Kerr black hole, 
which lies entirely out of the stable regime. 

FIG. 8. R{M/Mpf as a function of J/M^ and Q/M for (M, J) fluctuations. 
R is real and positive everywhere in the physical regime. R diverges at both 
limits of stability. 
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Figure 2 
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Figure 3 
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Figure 4 
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Figure 5 
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Figure 6 
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Figure 7 
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